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1.0 INTRODUCTION 


Satellite-borne, multi -spectral sensors such as found on the ERTS 
have produced a large volume of imagery v/hich, if it is to be useful, 
must be transmitted to earth and analyzed. The large volume of this 
imagery has in fact created problems in how to best transmit, store and 
analyze the images. A useful alternative to this procedure is to perform 
on-board -processing of the imagery for purposes of redundancy reduction 
or feature extraction. A successful on-board processing system would 
greatly reduce the load on transmission and storage systems. 

j 

The most promising candidate for on-board processing of images is a 
coherent optical data processing system. Such a system can effectively 
handle the large information content of imagery at high speeds. Evaluation 
and verification of coherent optical data processing for this application 
■involves theoretical as well as experimental knowledge of the characteristics 
and limitations of the system. This program had the objective of establish- 
ing techniques and a knowledge base for performing this evaluation and 
verification of an optical data processing system designed to reduce 

I 

redundancy in picture transmission and to detect specific image features. 

Theoretically derived Fourier transform characteristics for simple 
but representative two-dimensional images serve as a basis for predicting 
expected features of actual target images. Typical targets to be considered 
might include linear patterns representative of rows of agricultural crops. 
Large area patterns of varying’ shapes could represent large scale terrain' 
features. Targets of simple geometry v/ere analyzed first. Further refine- 
ment such as inclusion of random dot patterns of varying sizes or certain 
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curvilinear patterns is, however, necessary for adequate modeling of 
actual Imagery. 

Fourier transformation and spatial filtering of coherent optical 
images is readily accomplished in theory and in the laboratory. The 
effect of various parameters such as optical aperture, incidence angles, 
•the transparency assumption, and the thin lens approximation on resolution 
and performance of the optical data processing system were predicted and 

i 

tested. An exhaustive parametric evaluation was not done, however, those 
parameters shown to be most critical to the system operation were 

' i • 1 

evaluated thoroughly, d 

i I 
1 1 

Thus, the program consisted of analytical and experimental evaluation 

I 

of Fourier and correlation plane features in an optical data processing 
system for the purpose of redundancy reduction and pattern recognition. 
Figure 1 shows pictorial ly the main features of the program. 

The following tasks were to be accomplished in fulfillment of the 
proposed optical data processing research effort. 

Theoretical analysis . • 

• Determine Fourier and correlation plane features to be expected 
from simple, but representative, target model s- 

• Determine expected parametric limitations as related to feature 
extraction in the Fourier and correlation planes. 

• Study briefly detector geometry requirements as related to 

feature geometries. -• 

• Investigate methods for combined analysis of Fourier and 
correlation plane information. 
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Figure 1; Flow Chart for Optical Data 
Processing Study 
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Experimental verification 

-• Using GFE target film, verify actual Fourier and correlation 
plane features for several target classifications. 

• Compare experimental and theoretical features for representative 
targets. 

• Implement and evaluate candidate methods for simultaneous analysis 

i 

of Fourier and; correlation plane features. 



2;0 THEORETICAL ANALYSIS 


The development of coherent light sources has made possible the 
processing of data using optical methods. That such processing methods 
as Fourier transformation, correlation, convolution and filtering can 
be achieved optically is not obvious and needs further justification. It' 
is the purpose of this section to present a concise theoretical development 
showing that these processing methods can be implemented optically and 
to point out practical limitations to laboratory implementation of 
theoretically valid results. 

No attempt .is made to be complete in the description of the properties 
of coherent light. Rather the emphasis wilt’ be on techniques important 
to the goals of this contract. Additionally, the presentation is 
intended to be a concise reference useful to those persons familiar to 
some extent with coherent light systems and the mathematical development 
of Fourier transforms. All the information is available in scattered 
references, but nowhere is it available in as concise and readable (hope- 
fully) a form as presented here. 

Since light is a form of electromagnetic energy, a fundamental 
starting point for its description is Maxwell's equations. Through a 
series of assumptions (most of which can be imposed on a real system) and 
approximations the development leads to the Helmholtz equation and 
then the Fourier transforming properties of lenses. The assumptions and 
approximations are summarized in table form along with the result for 
easy reference. 

Fc'urier transformation leads logically to correlation, convolution 
and filtering. Theoretical- development of these operations is given. 


PAGE IS 
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The problems associated with the use of practical systems for 
data inputs recording and filtering are discussed briefly along with 
parametric limitations. 

Finally a brief theoretical introduction to advanced analysis 
methods which look at the fine structure of the optical correlation 
.plane is given. 
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MAXWELL*S EQUATIONS, THE HELMHOLTZ EQUATION AND THE 
FOURIER TRANSFORMING PROPERTIES OF LENSES 


Beginning with Haxv/elTs equations for a source free tnediunv (Assumption- §1) 


Vx ^ “ii . Vk H - ^ • Cl) 


It is convenient to use the temporal Fourier transforms of B & D 
to clear the time derivatives, i.e. 

ysc - 2.rriv3C'i) j YK - -zirCy C^) 


where : 


and 
^/v) and 



are Fourier transform pairs. 


Next we "assume a non-dispersive (uniform) medium so the dependence 
on V can be dropped. (Assumption n) 

Further if the medium is isotropic and linear we have (Assumption §3) 

D ^ • I r ^ ^ « (3) 


Otherwise e and are second rank tensors 

From the differential equations in (2) we obtain the 


VECTOR 

V £ -h 4TX^k £ - 0 

c^) 

HELMHOLTZ 

^ <v 

A . /V 


EQUATIONS 

Y H + 4-lT^k H ^ 0 



For a three dimensional medium (such as an optical system) solution 
of the Vector Helmholtz Equations involves simultaneous solution of six 
scalar differential equations subject to a set of boundary conditions. 
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Solution of Vector Helmholtz Equations is impossible except for a 

few (Impractical) casesi e.g., plane wave incident on perfect thin 

1 . - 

conducting half-plane. 


Scalar Approximation - Assume the optical phenomenon of interest 
is adequately described by a scalar quantity. Note: the optical 

phenomenon of interest ^is normally the magnitude of the electric field 
vector expressed as a function of X, Y and Z coordinates and either time 
or frequency. ' . 

The scalar approximation results in the 

t 

i 

1 

SCALAR fr\ 

HELMHOLTZ V' 4 O 

EQUATION 

A ^ • V 

where: 4 ~ ^ 

and • W 

■ 4 is real; 4 is generally complex 

If we assume a monochromatic field with temporal frequency 
= V*. 3 ' (Assumption w4) 


then 4 can be represented by 

a. at J j 

where: • is a general complex spatial amplitude function. 

For a plane wave , has the form 

f y y -f-T| ?:) 




- e 


Z.U-L (i^kX 


C?) 


- e 
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then . ^ zrri(A^X -i-A.u 

4>Cx.^,^,i)= ^ e -i 

The plane wave form of ^ Ch'ij'^i'^') satisfies the scalar Helmholtz 
equation if we impose the condition:. 

k - v/avenumber 

\\ \ =7 optical v/avelength. 

ii 

It can be shown that; a sum of planewaves also satisfies the scalar 

•tire, • 


Helmholtz equation, 

A ii 




where : 


and 


4'Cl‘i,^) ^ Z'/L(C 

wavenumber for nth planewave 


IOk ‘ — 

••v 


for all n 


. By taking this line of reasoning one step further, a continuous 
sum of planewaves also satisfies the scalar Helmholtz equation, i.e., 

ii ' 'L'n e, • A 




0 °) 


Ln) 


0^) 


0-3) 


Significance - Any optical function which can be represented as a 
sum of planewaves will satisfy the scalar Helmholtz equation. 


Assumptions to this point 



^/ViepitzMj) 


(ofnCAc Ff£UP) 

1. 

source-free 

1. 

adequately represented as a scalar 




quantity 

H. 

isotropic 


■ * . ' f 



2. 

monochromatic 

3. 

line ir 





3. 

can be represented as a summation 

4. 

uniform (non-dispersive) 


of plane-waves . 
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Propagation of Optical Fields 

Begin with the scalar Helmholtz equation given by- equation (6) 


v^cj> rl- — o . 

> . \ \ “Z-tr t- y^e- 

Define cf> == 3 

to separate temporal and .spatial dependence (This form also implies a 
monochromatic field) with being the complex amplitude of 

the field. ; 

Can verify that ^ also satisfies the scalar Helmholtz 

equation, i.e., ' ; • 

i 

* * , / 

Solve the boundary value problem for the field 
given C^i it o) as one of the boundary conditions.- The result is 
the Rayl ei gh“Sommerf el d integral 

where : ("I “ + ^""7 ^ 

Two approximations are now made to simplify (17) (See Figure 2) 


0‘i-') 

[is] 




Approximation 1 : ( ^ l) i.e.* distance s is much greater than the 
optical wavelength. 

Then we have the result known as Huygen's Principle 


~ co3[i,s) e 


ZiTt ks 




n 



that is, the field at (X, Y, 2) is a weighted sum of spherical waves 
emanating from 

Approximation 2 : (Fresnel approximation) at 2 = d (Figure 2), assume 

the size of the regions of interest in either plane is small compared to 
ds i.e., 

then: <i - 2 - 5 

CosCcijS) i 

J- ^ A 
5 ^ ? 

and the exponential term becomes approximately 

- s= dr I •)- 

^ 4- 4 . lu^Uyi -c^v>vwj 

oi-d 



12 


Since ^ sre small compared to d we can neglect the higher 

order terms and^ 


for exponent 



rx-^r-f N-i) " 

a-d 



Substituting (20) into (18) gives 




Equation (21) is the result of the Fresnel approximation and holds in 
the region close to the optical axis. 

Recognize that (21) is (See Appendix A) times the - 

Fresnel transform of VCx,if^o). 

The constant phase term is usually ignored and the result is that if the 
optical field is known in a given plane, it can be found in another plane 
by Fresnel transformation . 


COMMENT : The Fresnel approximation holds much better than expected. 

The reason is that for images with low spatial frequencies 

there is no longer a need to restrict the regions of interest to be close 

to the optical axis. Most images of interest will fall in this category . 

Fourier Properties of Lenses 

Transparency Assumption - Assume a thin lens has t(X, Y) transmission 
function consisting of the product of a pupil function and a pure phase 
term, i.e., 

= pU,^] 
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then the optical fields immediately before and after the lens are 
rel ated by ; 


TRANSPARENCY 

ASSUMPTION 






- Now the phase term in (22) is given by 
<^LXiv) - ^Trk(r>x-i') cUXi^)' 

where: n = refractive index of lens material 

k = wavenumber = 
define: focal length 

f i, 

then V 

A A-P 

If. we combine the constant term in (29) into the pupil function 
given in (22) we have 


c^^) 

(y-s) 

(y‘i) 
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(30) 


Equation (30) 1s the transmission function for a ”thin'^ lens . 

Now the focal plane response is found by 

Incident field = <xC^) at front side of lens 

After lens = 

t ♦ 

-ic. -In focal plane ' . 

Using given in (30) and after cancellation of some terms the 

result is; 


FOURIER ^ f- 1 

TRANSFORM- -JTT^ j ^ 




pi%) ^(f) 


(3i) 




This is the desired result. Except for the phase term outside the 
integral > equation (32) shows that the focal plane response is the 
Fourier transform of the product of the pupil function v/ith the incident 
field, i.e.. 
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COMPLEX FILTERS, .CONVOLUTION AND CORRELATION . 


COMPLEX FILTER GENERATION 



flGUf^E 3 V CorAPi^ 


Amplitude 1n the filter plane (See Figure 3) is the sum of the reference 
beam and the Fourier transform of - At Cx>^) » i.e., 

where; dtCXiy) includes any pupil function; 

■' to.,)- 




( 35 -) 


The photographic film in the filter plane is an energy sensitive medium 
thus the transmission of the developed negative is proportional to 
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fCif.j) = 1- 


Substituting in for 



we have 


1-^7 




- 1 — G^o 


A 

do* 





Through proper development of a positive transparency a transmission 
function proportional to | can be obtained. Either positive 

or negative transparencies will provide correlation results. 


COMPLEX FILTERING 




OUTPUT 

FuiA/p 





CoMF Fl LTF^jet /[/ G- 
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Figure 4 shows an arrangement for observing the cross correlation 

between images (Xi C^l‘^')and Ch. (Xiy) . It also gives the convolution 

of the two images. • . . - 

! 

■ The light amplitude inmedlately to the right of the spatial filter 
is given by: 




I 1 


^ yf) Tf)] if) 

iif [ A ^ e"'"] A»,*) 


•f". 



•After transforming by lens Z 2 . s equation (38) gives the follov/ing 
amplitude distribution in the output plane: 


where: 

THa/5 


r- 


+ i") ^ f(x y.\f) ^ 


ofiDsn 


I — 

k 


K= 

CMx>S^ /::./i-77oAj 


60 n r \ ^ / Cl co/e^/^A/i-77oAj 

—J 
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'■ Equation (40) is the desired result and for O.iCxty') = 
the resulting autocorrelation shows a strong peak at X" = ^o-PA » 

Because the spatial frequency content of most images is low, the resulting 
autocorrelation function appears simply as a bright spot. Energy 
detection for this so-called recognition spot is sufficient for distinguish- 
ing gross image differences. However, when comparing images that are 
quite similar, it is ne’cessary to examine the fine structure of the 
correlation function to aid in distinguishing betv/een these similar 
images. Examination of the fine structure of the correlation plane requires 
magnification, sampling and additional processing (probably by computer). 

An alternate approach to the vander Lug t method for generation of 
optical correlation functions described above is discussed below. 






Figure 5: Power Spectrum Recording 

The effect of the lens on the input image s{ir) is to produce the 
Fourier transform Photographic film in the Fourier plane 

responds to intensity and after development has a transmittance proportional 


to 
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Or t(^) 1s proportional to the power spectrum of the original image. 
If this pov/er spectrum transparency is viewed as shown in Figure 6 


AUTO COA^eiAT/OAJ 



Figure 6 : Autocorrelation 




Figure 7: Generation of Cross Pov;er Spectrum 
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Input 

SC%)~ Si X %) 4- ¥r %) 


Field in Fourier plane 




IX’f 

Developed film transmittance 





X ^ 

' XXf. 


_ rt-lrlui 


4— 


km e 


J 


; zirXj ^ 


~^>4 


? C^l - ! I ( vf) I 



Figure 8, Cross Correlation 

rl-f I ^ ^ 





•f- i?i^(‘2^) ^ 

f- /?z, ^ 


The relative strength' of the cross correlation is strongly dependent 
on the degree of angular alignment in making the cross-spectrum transparency 
in Figure 7. 
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PRACTICAL CONSIDERATIONS 

A number of parametric anal^'ses can be derived from limitations 
imposed by practical systems and from approximations made in the theo- 
retical development. Figures 9,. throughH are described below and provide 
parametric variations affecting the Fourier transforming properties of 
lenses . These figures were also included in Monthly Progress Report No. 5. 

Figure 9^ In order, for no vignetting to occur certain restrictions 
must be placed on the input, output, and lens .apertures. These relation- 
ships are depicted graphically where (max)/f is the (normalized) output 

1 

aperture and r(max)/f is the (normalized) input aperture. Given a lens 
with a particular F-stop, r and y? must li.e;.'.below that curve for no 
vignetting. 
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Figure 10 In the transform plane the spatial frequencies do not 
Increase linearly with the distance from the -optical axis. The per 
cent error in the linear' frequency approximation is shown as a function 
of the output aperture size (normalized). . 

Figure 11 and Figure 12 . These curves v;ere derived using the exact 
expressions for the obliquity factor and the distance from- the object 

I 

point to the image pointj (i.e., the Fresnel approximations were not 
made). Figure 11 shows the restrictions on input and output apertures 
(normalized) required to achieve a given maximum amplitude error in the 
transform plane. Figure 12 shows the phase term which multiplies the 
Fourier Transform expression for various values of focal length f. 

t S. 

Figure 1.3 Axial misalignment of the input plane from the lens 
results in a- phase error in the transform plane, shown for various 
values jd of maximum alignment error. 

• Figure 14 This figure shows the limits on frequency resolution 
for a perfectly aligned system. The limit is determined by the size of 
the aperture in the input plane. 

The following parametric considerations apply to complex filtering . 

Minimum reference beam angle - With reference to Figures 3 and 4 
and equation (40), if both and have width, W, 

then separation of terms in the- output plane of Figure 3 will be 
achieved if 

■ if ( 4 ,) 
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Proof: In equation (40) we have 

TOO ZERO-ORDER TERMS (centered at X, Y = 0, 0} 

I 

1st term has width W 
2nd term has width 3W 

CONVOLUTION TERM (centered at X, Y = , 0 } 

has width 2 W 

CORRELATION TERM (centered at X, Y = j 0 ) 

has width 2W 
Separation achieved if 

^„-PA-u/> iw 

or 

A = 0 > f w 

and finally with reference to Figure 4, a minimum aperture for lens I2 
can be established to ensure that the correlation and convolution terms 
of equation (40) fall within the lens aperture. 

Film Resolution Requirements - For interference of two plane waves 
with angular separation oC , the interference fringes are separated by 

d - _A_ (-4^^ 

s I r> c< 

Figure 15 ■ shows a plot of the reciprocal of this equation for 
62-zg /9 . It shows how resolution in cycles per millimeter 

Increases as the reference angle increases. 



’.. j {^ereie^/oae' dET-if^ /h/\/(SC£ 

* . . r * ♦ ‘ 

',^&i{/e£_ /p , . . ^/i. A^i , oca r/ oa/ ^ o<k/^ /. //ti iV 
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Effect of inherent phase term in Fourier transform expression - 
Even 1f the aperture restrictions are such that the amplitude error in 
the Fourier transform expression may be neglected and the linear 
frequency approximation is validj there still exists a multiplicative 
phase term given by 

Wan 


£ 


!i 


which may not be neglected. However, the following development shows 
that this term is of no consequence in the operation of the optical 
correlator. [It does however affect the coherent imaging process.] 

Only the last term in the expression for T 'will be considered 
since it is the terra which will result in the correlation term in the 
output. To include the multiplicative phase factor this term should be 
written: 


-JL 

iX-f 






C 




When this is multiplied by the Fourier transform expression 






a) 

Af-i 


C 


(4s) 


( 44 ) 


( 47 ) 


these two additional phase terms cancel resulting in the original .. ■- 
expression. Since a longitudinal positioning error of the input trans- 
parency results in a similar phase term, no error will result, providing 
the filter is recorded and all input signals are placed in the same 
longitudinally displaced position. 

A similar statement may be made concerning lateral displacement of 
the inpiic image. Since this lateral displacement results in the Fourier 


original?^® 
OF POOR eOALm 
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transform expression being multiplied by a phase factor 



Izrrfd 


where d is the lateral displacement and is radial distance in the 
transform plane, then aslong as the recording and correlation processes 
take place using the same laterally displaced position, no error results. 

This is true because tHe phase term times its .complex .conjugate equals 
unity. 

No analagous general statement can be made concerning angula.r 
positioning and rotation since these involve other than phase variations. 

Two other types of displacements may occur in the input plane. These 
are angular displacements and rotations. Angular displacements result 
in theimage plane no longer being perpendicular to the optical axis; 
whereas, rotations are defined to be angular movement about the optical 
axis v/hich maintains perpendicularity. In general, the effects of these 
displacements are strongly dependent on the particular image function. 

As an example, the periodic signal cos^zTryUxx) after being 
angularly displaced by the angle 0 becomes the new signal 

cos ^e) I^b) 

In effect the "lines" of the signal have been "compressed" resulting 
in an expansion in the transform plane. 

Effects of rotations in the input plane depend strongly on the signal 
since those signals with circular synmetry being correlated with themselves 
are immune to rotations while those with straight-line features may be 
affected greatly. Effects of rotation may be found for particular signals 


ORIGINAL PAGE IS 
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g(x, 3 ') and h(x,y) but little can be said in general since rotations may 
either increase or decrease the degree of correlation between two 
unlike signals. It might be noted that the effects of rotations '(uni ike 
lateral V longitudinal, and angular displacements) are the same (for a 
given rotation $) whether the rotation takes place in the input or 
filter plane. 

I 

Filter displacements are much more critical than are image displace- 
ments between recording and correlation. In fact the above parametric 

i 

. ' t 

effects of small image displacements assume that the filter is perfectly 

1 i 

aligned. So critical is this alignment that it normally is accomplished 

I 

1 

• with the aid of; a microscope. The following development treats 

I 

the effects of filter misalignments and is derived from A. Vender Lugt 
(Applied Optics, V.6, N.7, duly 1967). 


Figure -16 . Shows relative correlation intensity for small lateral 
displacements of the filter for three ratios of image maximum size to 
focal length. 

Figure d7. Shows relative correlation intensity for small longitudinal 

\ 

displacements of the optical filter as a function of image max. size to 
focal length ratio for three nominal displacements. 


To determine the effect on performance of angular displacement of 
the filter, some previous results may be used- For a tilt of the 
longitudinal displacement of the filter is given by 
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v/here y^Onax) 1s the maximum extent of the Fourier Transform plane. 
Using this^Z in the longitudinal displacement equation yields 
performance for different Q {See Figure 18). 



3.0 EXPERIMENTAL ANALYSIS 


Experimental work on the program was impeded on several occasions 
because of construction and modification to the electro-optics laboratory. 
However, the major experimental tasks were accomplished. Improved 
results could have been obtained with more sophisticated auxiliary 
optical equipment such as precision positioning and alignment fixtures. . 

From the aerial photographs supplied by NASA, segments in 35 mm. 
format were selected which represent different types of terrain. For 
example, from the supplied photographs, several textures of cultivated 
land, orchards in various stages of development, forests, housing, motor 
vehicles and combined water-land areas wer^:'pbtained. Since multiple 
shots of all photos were provided, the procedure' was to simply cut example 
areas from one copy of each different photo. This method is simpler 
than a photographic reduction process and further has the advantage of 
retaining the resolution qualities of the larger format film. The only 
problem is that because of larger images the spatial frequency content 
will be more concentrated near the optical axis In the optical transform 
plane. 

A total of eighteen representative segments were selected and 
optically transformed to check for significant differences and similarities 
among the various spectra. Photographs of the transforms were taken and 
enlarged/printed to facilitate comparison. Additionally, autocorrelation 
functions were generated optically and photographed for each of the 
candidate images. The method used for generating the correlation functions 
was to transform optically a photographic transparency of the power spectrum 
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for each image. This method eliminates alignment problems encountered 
in the Vender Lugt method for generating autocorrelation functions. 

i 

Cross-correlation functions may be generated in. a similar fashion; 
however, there will be an alignment problem during recording of the 
cross-power spectrum. 

A number of correlation experiments using the Vander Lugt method 
were conducted with the result being essentially the same for all 
experiments. That is, no recognition spot at all was obtained with any 
image except the one used to make the spatial filter. Even the auto- 
correlation spot was very weak -in some cases, an effect most likely. 

attributable to the previously mentioned alignment difficulties. The 

/ 

T '-.4 - 

cross-correlation failed even for segments'cut from the same larger 
image format (eg., two segments of a peach orchard). Enlarged photographs 
of selected aerial imagery, power spectra and autocorrelation functions 
are shown on the following pages. The need for more accurate alignment 
control is indicated by the results. 
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Figure 19. Peach Orchard (0132055, #3) a) image b) power spectrum 
c) autocorrelation 
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Figure 20. Peach Orchard (0132055, #1) a) image b) power spectrum 
c) autocorrelation 
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Figure 21. Plowed field (01320063, #1) a) image b) power spectrum 
c) autocorrelation 
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Figure 22. Pasture (01320063, #3) a) image b) power spectrum 
c) autocorrelation 
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Figure 23. Orchard (0132055, #2) a) image b) power spectrum 
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Figure 24. Peach Orchard (01320120, #2) a) image b] power spectrum 








Figure 25. 


Plowed field (01320120, #4) 


a) image b) power spectrum 
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Figure 26. Sparsely populated orchard (01320120, #1) a) image b) power 
spectrum 
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Figure 28. Peach orchard, various tree sizes {01320063, #2) a) image 
b) power spectrum 
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4.0 COMPUTER ANALYSIS 


Coinputer programs -for computing the two-dimensional Fourier transform 
and providing a three dimensional perspective plot were generated. The 
purpose was to Investigate the feasibility of generating computer results 
for comparison with expected optical transforms. The programs .were 
successfully operated, the major limitation being, as expected, limited 
resolution because of core and time requirements. However for relatively 
low resolution Imagery the program does provide transform and correlation 
results. In order to use the programs effectively in analysis of images 
it is necessary that the Imagery be scanned and digitized. Equipment 
for performing the scan and digitization as .well as transform and filtering 
operations is available commercially at significant cost. 

For this project, fortran programs were generated using the standard 
IBM FFT subroutine HARM which accepts two dimensional or even three 
dimensional data sets. Additionally a program suitable for operation on 
smaller computers v/hich utilizes manipulation of two dimensional data in 
a one dimensional FFT algorithm was implemented. The procedure is to 
first transform the rows of a two dimensional data matrix and then 
transform the columns. 

Arrays up to 128 x 128 were transformed and plotted on the IBM 370 
and Calcomp Plotter. CPU time and core requirements for the array size's 
were approximately as follows: 

array size CPU .time core 

128 X 128 30-40 sec 394K 

64 X 64 • , 25-35 sec 150K 
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Example outputs are shown below: 

. Figure 29. Calculation and plot of magnitude of Fourier Transform 
of square aperture (64 x 64 array). 

Figure 30. Calculation and plot of magnitude of Fourier Transform 
of circular aperture (128 x 128 array). Larger size array is used for 

I 

more accuracy in image definition; however, only the central (64 x 64) 

portion of the transform is plotted. 

Figure 31. Expanded quadrant plot of the spectrum for a circular 

aperture. , j 

M •! 

Figure 32. Spectrum for a narrow slit parallel to the y axis. 
Figure 33. Spectrum for a narrow slit parallel to the x axis. 
Figure 34 and Figure 35. Optically generated spectra for the 
rectangular and circular apertures. 

The computer program listings are given in Appendix B. 
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5.0 CONCLUSIONS 


The theoretical development for optical processing is based on a 
number of assumptions and approximations. Examination of the parametric 
effects of these assumptions and approximations indicates very little , 
error in the presentation of transforms and correlations. However » other 
parameters, specifically errors in alignment, can be shov/n to have a 

I 

strong effect on the success of optical processing techniques. This 


conclusion was predicted from parametric considerations and verified 
experimentally. 

A major problem in the experimental implementation of optical 
correlators, in addition to the alignment problem, is due to the fact 
that the spatial frequency content of most imagery is clustered near 
the optical axis giving a very bright central spot. Any optical detection 
medium such as film then has the requirement of responding to optical 
patterns v/ith a wide dynamic range. Saturation effects are almost a 
certainty if- exposure is increased to detect the weaker pattern regions. 

Further efforts at application of coherent optical processing 
techniques in other than a laboratory environment requires solution of 
practical problems. 
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Figure 34. 


Figure 35. 
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Optically generated spectrum for a rectangular aperture. 



Optically generated spectrum 


for a circular aperture. 
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APPENDIX A:' FRESNEL TRANSFORMS 


A Fresnel function is defined as 

+ (x):i= ^ - e 




where ^ is positive and real. Its Fourier transform is found from 

'1 X j 

F(a] e 

I * j 

Equation (A2) can bej evaluated using Fresnel integrals to obtain 


F(/&) = 


cn^ya. 






Similarly we have the Fourier transform pair 




-tTT-CL'E 


r-.l, 

<^=^> ~-7T ^ 


From the convolution theorem and the above 


€ ¥ e a 




Now define the Fresnel transform of g(x) 






V 

where g(x) is the Fresnel transform of g(x). Using property (A5) it is 


easy to show that 






(A6) and (A7) define a Fresnel transform pair. 
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Comparison of Fourier and Fresnel transforms 


If we Fourier transform both g{x) and its Fresnel transform g{x) 
v/a get ■ • 


I _ 




m 


gC%) 

m 


Substituting (A6) Into (A9) gives 


t(J) =1 




'i' e 






nteo 


or 


& (fi) = e 


(All) 


Equation (All) says that a Fresnel transform in the spatial domain 
corresponds to a multiplicative phase -factor in the spatial frequency 
domain. 

'Fresnel transforms are useful in describing propagating optical 
fields. and in the study of Fresnel diffraction. 



APPENDIX B: 
COMPUTER LISTINGS 
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C PROGRAM USIMG SUBROUTINES FFT AMD PL0T3D 

COMPLEX A (64 5 ,.U, W , T r F { 6^ , 64 ) t C ( 64 , 64) 

DIMENSION OUTBUFI 64) ,MASK( 2000) , VERTEX! 16) 

N=64 U 

M=6 

C(I?J) IS INEtlALLY THE COMPLEX N X M ARRAY TO BE . 
FOURIER TRANSFORMED 

Ad) IS THE ONE DCMENSCOMAL N-ELEMENT ARRAY ON WHICH 
SUBROUTINE FFT OPERATES 

UtWjAND T ARE VARIABLES OF SUBROUTINE FFT 
OUTBUF,MASK» AND VERTEX ARE ARRAYS OF SUBROUTINE PLOT3D 
FORM INPUT ARRAY FOR SQUARE APERTURE 
DO I I = 1*M ,, 

DO 1 J=1,N 

Cd,J)=0.0 i ■ ■■ 

IF ( JoGE,29o AND. JoLE* 36.ANDd oGE.29»ANDd .LE.36 ) 
«CCI»J)=L.O 

1 CONTINUE 

PERFORM ONE DIMENSIONAL FFT ON ROWS OF C{I,J> AND 
STORE RESULT IN Fd,J) 

00 4 1=1, N 
00 2 J=1,N 

2 A(J)=Cd7J) 

CALL FFT(A,M,N) 

DO 3 K=1,M 
Kl=N-K+i 

3 F( I,K1 )=A(K) 

4 CONTINUE 

EACH OF THE N-ELEMENT ONE DIMENSIONAL TRANSFORMS NOW 
FORMS A ROW OF F(I,J) 

TAKE TRANSFORM OF COLUMNS OF F(I,J) AND STORE RESULT 
IN CdvJ) 

DO 5 J=l,N 
DO 6 1=1, N 

6 A( I )=Ft I, J) ■ 

CALL FFT(A,H,N) 

DO 7 K=1,N 
K1=N-K+1 

7 C(Kl,J)=A{K) 

5 CONTINUE 

TRANSFORM IS NOW STORED IN C(I,J). TAKE ABSOLUTE VALUE 
OF EACH COMPLEX ENTRY 
DO 8 1=1, N 
DO 8 J=1,N 

8 C( I, J)= CABS!C{ I, J) ) 

C REARRANGE FOR PROPER OUTPUT AND STORE IN ARRAY F(I,J) 

ND2=N/2 
ND2Ml=ND2-i 
N02P1=ND2+1 
ND2P2=ND2+2 
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' DO 71 I = ltND2Pl 
DO 71 J=1,ND2P1 

71 F ( H-ND2M 1, J+N02M1 }=C( f,J} 

DO 72 I=ND2P2,iM 

DO 72 J=1,ND2P1 

72 F(I-ND2P1,J+ND2M1)=C{I,J} 

DO 73 I = .\D2P2,N 

DO 73 J = r.JD2P2,N 

73 F(I-ND2PI,J'*-ND2P1) = C{ I»J) 

DO 74 I=lrMD2Pl 

DO 74 J=N02P2tN 

74 F( H-ND2M1, J-ND2P1 )=C( I, J ) 

C NORMALIZE VALUES IN F(I,J) TO VALUE OF ONE 

FMAX=0.0 . . • 

DO 21 1=1, N i 

DO 21 J=1,N 

IF {FMAX.GEaCA8S(F(I, jn } GO TO 21 
FMAX=F(I,J) 

21 CONTINUE i 

DO 22 1=1, N ; 

DO 22 J=l,N 

F( I, J)=F( I, J)/FHAX 

22 CONTINUE • ‘ 

C PLOT THE ARRAY f!{I,J) USING SUBROUTINE PLOT3D 
00 20 NLINE=1,N , 

DO 10 NP0INT=1,N 
OUTBUFINPOINT)=F{NPOINT,NLINE) 

10 CONTINUE ' 

CALL PL0T3DI 10 10,0 oO,OUTBUF, 0.0, 0.075 ,4. 0,-0. 075, 
INLINE, 64, -45.0, -30.0, 5. 0,3.0, 10. 0 , MASK , VER TE X ) 

20 CONTINUE 

C DRAW FRAME AROUND PLOT USING SUBROUTINE FRAMER 
CALL FRAMERO, VERTEX, MASK) 

C SIGNAL THAT PLOT IS COMPLETE 
CALL PLOnO. 0,0.0, 999) 

STOP 

END 
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SUBROUTINE FFTCAsM^N) 

COMPLEX A(X'),U»V.ST 

NV2=M/2 . 

. • •; 

DO 7 I=lrNMI 
IFU.GE.J) GO TO 5 
T=--A{J) 

A(J)=A(I) 

AU) =T 

5 K=NV2 ir 

6 IF(K.GE.J) GO TO-7 
J = J-K 

K=K/2 
GO TO 6 

7 J=J-t-K 
Pl=3. 14159 
DO 20 L=li>H 
LE=2=i^’{'L 
•LEl=LE/2 

u=(i«o,o.) 

W=CMPLX(COS (PI/LEl} , SIN(PI/LE1 ) ) 
00 20 J=lfLEl 
DO 10 I=J,N»LE 
IP=I+LEi 
T=A( IP)«U 
Af IP )=A( I )-T 
10 At [)=A{I ) -f-T 
20 U=U«W 
. RETURN 
END 
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PR0G9A?-'. USING SUBROUTINES HARM AND PL0T3D , 
DIMENSION A{3E768) ,M{3) , I iNV ( 4 096 ) » Y ( 12 B p 1 2B-) r 
’!'Xil2 8,120)TS{ 4 096 ) i I VAL( 128) , MASK (2000 } » VcRTEXt 16 ) 
«0UTBUF(128) 


s 


N=128 . - ■ 

MSQD=N'-:"^s2 
NSQDT2= 2>!=NSQD 

FORM INPUT ARRAY X(lrJ) TO BE FOURIER TRANSFORMED 
AND COLUMNiZe TO ARRAY All) FOR INPUT TO SUBROUTINE 
HARM . - • 


K = 1 

DO 1 J=ltN 
DO 1 1=1 tN 

FORM INPUT ARRAY FOR CIRCULAR APERTURE 


X( I,J}=0«0 

. IF (JcGE«62.AND.J.LE.6 7i.AND.I.GE«57 


❖ *AND*I *LE .72) iX(l7J) = U0 
IF ( ( J.EQ.57.0R. J. EQo 72) . AND . I . GE . 62 » AND . 
»!'I<.LE.67) X(I,J) = 1.0 

IF (UoEQ.58.0R. J.EQ.71) . AND . K GE o 60 . AND » 
>}=KLE«69) x(rtJ:) = i.o 

IF ( (JoEQ.59.0R.J.EQ,70) . AND. I . GE»59, AND. 
=MoLE.70) X{ItJ') = 1.0 
IF t ( J.EQ.60.0'R, J.E0.61.0R,J.EQo68.0Rc 
=?'JoEQ.69) . AiND.Il .GE. 58. AND. I .LE.71 ) X(I»J)=loO 
A(K)= X( I , J) I 
A{K+l)=0o0 
K=K+2 


1 CONTINUE . 

C SET UP PARAMETERS FOR SUBROUTINE HARM 

IFSET=1 
H(l)=7 
M(2)=7 

. M(3)=0 . ' . • 

CALL HARMtAjMt INV, St IFSET 7 IFERR) 

C CHECK FOR ERRORS 

IF (IFERR.NE.O) GO TO 35 
■ GO TO 60 

35 WRITE(6,50) IFERR 

50 FORMAT ( * 1 * , 30|( / ) t 56X , » IFERR= ' , I 2 ) 

STOP 

C . PUT THE COMPLEX TRANSFORMED ARRAY A(l) INTO NEW 
C ■ MAGNITUDE ARRAY A(I) 

60 DO 65 I=1,NSQDT2,2 

65 A{(l + l)/2)= SQRT( (A( I) )^*2+<A(I + l) 

C PUT ARRAY Ad) INTO TWO DIMENSIONAL ARRAY Y(I,J) 

1=0 • ' - 
DO 70 L=UN 
DO 70 J=1,N 
1 = 1 + 1 

70 Y{J,L)= A( I) ' 

C REARRANGE FOR PROPER OUTPUT AND STORE IN ARRAY X(IiJ) 

ND2=N/2 
ND2Hl=ND2-l 
ND2Pl=MD2+l 
ND2P2=ND2+2 
' DO 71 I=ltN02Pl 
DO 71 J=l,NU2Pi 
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71 xn+ND2MlrJ+ND2Ml)=Y(I,J] 

DO 72 I = r-}D2?2tn 

DO 72 J=1,HD2P1 

72 X( I-ND2pit J+ND2H1 )=YtI ,J) 

DO 73 I = tN!D2P2,N 

DO 73 J=ND2P2,N 

73 X{ I-ND2PI, J--Nn2Pl)=Y( I rJ) 

DO 74 I=1,ND2PI 

DO 74 J=ND2P2rN 
.. 74 Xt I+ND2MI,J-:JD2P1)=Y(1 ,J) 

C NORMALIZE VALUES IN X(1,J) TO VALUE OF ONE 

XMAX=0.0 
DO 21 1=1, N 
DO 21 J=1,N 

IF (XMAXo GE.XC I , J ) ) GO TO 21 
XMAX=X( I , J) 

21 CONTINUE 

DO 22 I=1,N 
DO 22 J=l,N 
X(I, J)=X{ I , J)/XMAX 

22 CONTINUE 

PLOT THE CENTRAL 64 X 6A ELEMENTS OF THE 128 X 128 
ARRAY XI I, J) USING SUBROUTINE PLOT3D 
DO 20 1=33,96 
DO 10 J=33,96 
NLINE=I-32 
NPOlNT=J-32 
OUTBUF(NPOINT)=X( J, I ) 

10 CONTINUE 

CALL PLOT 3D { 1010,0, OjOUTBUF, 0.0,0.075 ,4.0,-0.075, 
=«‘NLINE, 64, -45.0, -30.0, 5. 0,3.0 , 10,0 , MASK, VERTEX ) 

20 CONTINUE 

CALL FRAHER{3, VERTEX, MASK) 

CALL PLOTtO, 0,0. 0,999) f 
STOP 
END 


OP.POOK pUALIIY 




ELECTRICAL 


ENGINEERING EXPERIMENT STATION 

AUBURN UNIVERSITY 

AUBURN, ALABAMA 



